Mon. Not. R. Astron. Soc. 000, 000-000 (0000) Printed 5 February 2008 (MN WF$t style file v2.2) 



High eccentricity planets from the Anglo- Australian 
Planet Search 



Or 



ABSTRACT 



Hugh R. A. Jones 1 , R. Paul Butler 2 , C.G. Tinney 3 , Geoffrey W. Marcy 4 , Brad D. 
^ ; Carter 5 , Alan J. Penny 6 ' 7 , Chris McCarthy 8 , Jeremy Bailey 9 

^— ^ 1 Centre for Astrophysics Research, University of Hertfordshire, College Lane, Hatfield, Herts ALIO 9AB, UK 

Carnegie Institution of Washington, Department of Terrestrial Magnetism, 5241 Broad Branch Rd NW, Washington, DC 20015-1305, USA 
J.^ ' 3 Anglo- Australian Observatory, PO Box 296, Epping. 1110, Australia 

4 Department of Astronomy, University of California, Berkeley, CA, 94720, USA 
^5 1 5 Faculty of Sciences, University of Southern Queensland, Toowoomba, QLD 4350, Australia 
6 Rutherford Appleton Laboratory, Chilton, Didcot, Oxon OX 11 0QX, UK 
7 SETI Institute, 515 N. Whisman Road, Mountain View, CA 94043, USA 

8 Department of Physics and Astronomy, San Francisco State University, San Francisco, CA 94132, USA 
9 Australian Centre for A strobiology, Macquarie University, Sydney NSW 2109, Australia 

> 

in 

co 

CO 
CO 

o 

\^ | We report Doppler measurements of the stars HD 187085 and HD 20782 which indi- 

, cate two high eccentricity low-mass companions to the stars. We find HD 187085 has a 

Jupiter-mass companion with a ^1000 d orbit. Our formal 'best fit' solution suggests 
an eccentricity of 0.47, however, it does not sample the periastron passage of the com- 
panion and we find that orbital solutions with eccentricities between 0.1 and 0.8 give 
only slightly poorer fits (based on RMS and xl) an d are thus plausible. Observations 
^—j | made during periastron passage in 2007 June should allow for the reliable determi- 

nation of the orbital eccentricity for the companion to HD 187085. Our dataset for 
HD 20782 docs sample periastron and so the orbit for its companion can be more reli- 
ably determined. We find the companion to HD 20782 has M sin i = 1.77±0.22 M Jup , 
an orbital period of 595.86±0.03 d and an orbit with an eccentricity of 0.92±0.03. The 
r% \ detection of such high-eccentricity (and relatively low velocity amplitude) exoplanets 

^ . appears to be facilitated by the long-term precision of the Anglo- Australian Planet 

Search. Looking at exoplanet detections as a whole, we find that those with higher 
eccentricity seem to have relatively higher velocity amplitudes indicating higher mass 
planets and/or an observational bias against the detection of high eccentricity systems. 

Key words: planetary systems - stars: individual (HD187085) (HD20782), brown 
dwarfs 



1 INTRODUCTION 

The Anglo-Australian Planet Search (AAPS) is a long-term 
radial velocity project engaged in the detection and mea- 
surement of extra-solar planets (hereafter, shortened to ex- 
oplanets) at the highest possible precisions. Together with 
programmes using similar techniques on the Lick 3 m and 
Keck I 10 m telescopes (Fischer et al. 2001; Vogt et al. 2000), 
it provides all-sky planet search coverage for inactive F, G, 
K and M dwarfs down to a magnitude limit of V=7.5. So far 
the AAPS has has published data for 23 exoplanets (Butler 
et al. 2001; Carter et al. 2003; Jones et al. 2002a,b, 2003; 
McCarthy et al. 2005; Tinney et al. 2001, 2002a, 2003, 2004, 
2005, 2006) as well as non-detections of planets claimed else- 
where (Butler et al. 2001; Butler et al. 2002a). 



The AAPS is carried out on the 3.9m Anglo- Australian 
Telescope (AAT) using the University College London 
Echelle Spectrograph (UCLES), operated in its 31 lines/mm 
mode together with an I2 absorption cell. UCLES now uses 
the AAO's EEV 2048x4096 13.5ttm pixel CCD, which pro- 
vides excellent quantum efficiency across the 500-620 nm I2 
absorption line region. From 1998 until 2002, the AAPS al- 
location was around 20 nights per year, from 2002 to 2004 
it increased to 32. In 2005 the AAPS has an allocation of 
64 nights per year. Despite this search taking place on a 
common-user telescope with frequent changes of instrument, 
we have historically achieved a 3 m s" 1 precision down to the 
V = 7.5 magnitude limit of the survey (e.g., fig. 1, Jones et 
al. 2002a), for stars with low 'jitter'. In addition to our pub- 
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Figure 1. UCLES spectrum of the region of the CaHK line for 
HD 20782 (middle, log R' HK = -4.93, G3V) and two comparison 
objects of similar spectral types and activities HD 17051 (top, 
log R' HK = -4.81, GOV) and HD 142 (bottom, log R' nK = -4.95, 
GOV). The spectra are taken from Tinney et al. (2002). 



lished errors, there are additional sources of error collectively 
termed 'jitter'. The magnitude of the 'jitter' is a function of 
the spectral type of the star observed. Wright (2005) gives 
a model which estimates the jitter for a star based upon 
its activity. Our continuing attainment of 3 m s -1 single- 
shot precision for suitably inactive stars is demonstrated in 
McCarthy et al. (2004) and Tinney et al. (2005). 

Our target sample, which we have observed since 1998, 
is given in Jones et al. (2002b). It includes 178 late (IV-V) F, 
G and K stars with declinations below ~ —20° and is com- 
plete to V<7.5. We also observe sub-samples of 20 metal-rich 
([Fe/H]>0.3) stars with V<9.5 and 7 M dwarfs with V<11 
and declinations below ~ —20°. The sample has been in- 
creased to more than 250 solar-type stars to be complete to 
a magnitude limit of V=8. Where age/activity information is 
available from log .Rj^indices (Henry et al. 1996; Tinney et 
al. 2002b; Jenkins et al. 2006) we require target stars to have 
log -Rjjk < ~4-5 corresponding to ages typically greater than 
3 Gyr. Stars with known stellar companions within 2 arcsec 
are removed from the observing list, as it is operationally 
difficult to get an uncontaminated spectrum of a star with 
a nearby companion. Spectroscopic binaries discovered dur- 
ing the programme have also been removed and are reported 
elsewhere (Blundell et al. 2006). Otherwise there is no bias 
against observing multiple stars. The programme is also not 
expected to have any bias against brown dwarf companions. 

The observing and data processing procedures follow 
those described by Butler et al. (1996, 2001) and Tinney et 
al. (2005). All data taken by the A APS to date have been 
reprocessed through our continuously upgraded anlaysis sys- 
tem. Here, we report results from the current version of our 
pipeline. 



2 STELLAR CHARACTERISTICS OF 
HD 187085 AND HD 20782 

The characteristics for HD 187085 and HD 20782 are sum- 
marised in Table Q Houck (1982) assigns spectral types of 
GOV and G3V to HD 187085 and HD 20782 which are con- 
sistent with information derived for and from the Hippar- 



Table 1. The stellar parameters for HD 187085 and HD 20782 
tabulated below are taken from a range of sources: spectral types 
from Houck (1982), activities from Henry et al. (1996) and Tinney 
et al. (2002); Hipparcos measurements from (ESA 1997); lumi- 
nosities, masses, temperatures, metallicities and space velocities 
from Nordstrom et al. (2004, abbreviated as N04 in the table, the 
errors given for Nordstrom et al. values are the average of disper- 
sions, no upper confidence limit is given for the age of HD 20782 
by Nordstrom et al.) and Fischer & Valenti (2005, abbreviated as 
FV05 in the table); jitter values are taken from Wright (2005). 



Parameter 

Spectral Type 
log R' HK 
Hipparcos N ^, s 
Hipparcos a 

log (Lstar/LQ) 

M star /M - N04 
M star /M - VF05 
T e fl (K)- N04 
T eff (K) - VF05 
[Fe/H] - N04 
[Fe/H] - VF05 
U,V,W - N04 (km s" 1 ) 
Age (Gyr) - N04 
Age (Gyr) - VF05 
jitter (m s^ 1 ) 



HD 187085 HD 20782 



GOV 
-4.93 
113 
0.0013 

1 q+0.55 

1 16 + S ' 06 
1 - 1D=C ().<)8 

6011±76 
6075±30 
0.05±0.10 
0.05±0.03 
11,-20,-14 
3.9±fl 
3.3±°;« 
2.80 



G3V 
-4.91 
191 
0.0009 
1.25±o. ? 

0.90±g;°* 
1.0±0.03 
5636±76 
5758±30 
-0.16±0.10 
-0.05±0.03 
-38,-61,-2 

13.0.4.5 
7 1+ 1 - 9 
4.11 



cos mission (ESA 1997). Both stars have been included in 
large-scale studies of nearby solar type stars. Nordstrom et 
al. (2004) included them in a magnitude-limited, kinemati- 
cally unbiased study of 16682 nearby F and G dwarf stars. 
Valenti & Fischer (2005) included them in a study of the 
stellar properties for 1040 F, G and K stars observed for the 
Anglo- Australian, Lick and Keck planet search programmes. 
Valenti & Fischer used high signal-to-noise echelle spectra 
originally taken for template radial velocity spectra and 
spectral synthesis to derive effective temperatures, surface 
gravities and metallicities whereas Nordstrom et al. (2004) 
used Stromgren photometry and the infrared flux method 
calibration of Alonso et al. (1996). Both studies use Hip- 
parcos parallaxes to convert luminosities in order to make 
comparisons with different theoretical isochrones to derive 
stellar parameters. To determine, stellar masses and ages 
Valenti & Fischer, use Yonsei-Yale isochrones (Demarque et 
al. 2004) and Nordstrom et al. (2004) use Padova isochrones 
(Giradi et al. 2000; Salasnich et al. 2000). Both sets of de- 
rived parameters agree to within errors. 

Both HD 187085 and HD 20782 have moderate activity 
indices log -Rjjk ~ ~4-9 from Henry et al. (1996). For the case 
of HD 187085 we also have a confirmation measurement from 
our own CaHK programme (Fig. [T]and Tinney et al. 2002) 
Furthermore there is no evidence for significant photomet- 
ric variability for either star in measurements made by the 
Hipparcos satellite. Combining Hipparcos astrometry with 
their radial velocities, Nordstrom et al. (2004) determine U, 
V, W space velocities. Based on Edvardsson (1993) these 
ages are consistent with the derived ages in Table 
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Table 2. Relative radial velocities are given for HD 187085. Ju- 
lian Dates (JD) are heliocentric. RVs are barycentric but have 
an arbitrary zero-point determined by the radial velocity of the 
template. 



JD RV error 



(-2451000) 


(m s 1 ) 


(m s 


120.9170 


-12.1 


5.3 


411.0753 


-2.5 


6.5 


683.1693 


16.1 


5.7 


743.0494 


9.3 


5.5 


767.0046 


8.8 


4.7 


769.0652 


5.3 


1.4 


770.1153 


5.3 


5.2 


855.9477 


8.6 


8.9 


1061.2140 


-4.5 


5.3 


1092.0512 


-28.2 


5.9 


1128.0267 


-25.9 


5.1 


1151.0146 


-15.4 


6.8 


1189.9230 


-6.4 


4.2 


1360.2816 


-7.3 


1.6 


1387.2168 


-8.7 


3.5 


1388.2355 


-12.1 


3.7 


1389.2664 


-12.9 


4.3 


1422.2121 


-6.9 


4.0 


1456.0917 


-7.9 


1.9 


1750.2587 


13.5 


3.8 


1752.2311 


11.9 


3.8 


1784.2071 


19.6 


3.2 


1857.1218 


-4.2 


3.7 


1861.0168 


9.2 


4.5 


1942.9842 


17.2 


4.2 


1946.9260 


17.9 


3.2 


2217.0666 


-12.1 


3.7 


2245.0429 


-24.0 


5.9 


2484.3002 


-3.0 


3.7 


2489.2600 


-5.2 


3.4 


2507.1966 


-0.7 


2.8 


2510.2154 


-2.1 


3.3 


2517.2619 


-3.5 


4.1 


2520.2747 


-10.4 


4.2 


2569.0854 


-5.8 


3.9 


2572.1667 


-8.3 


4.5 


2577.0269 


2.1 


3.0 


2627.9629 


9.8 


5.8 


2632.0744 


6.7 


3.4 


2665.9534 


18.9 


3.7 



3 ORBITAL SOLUTION FOR HD 187085 

The 40 epochs of Doppler velocity measurements for 
HD 187085, obtained between 1998 November and 2005 Oc- 
tober, are shown graphically in Fig. [3] and listed in Table 
121 Our velocity measurements are derived by breaking the 
spectra into several hundred 2 A chunks and deriving veloc- 
ities for each chunk. The velocity error, given in the third 
column and labelled 'error' is determined from the scatter 
of these chunks. This error includes the effects of photon- 
counting uncertainties, residual errors in the spectrograph 
PSF model, and variation in the underlying spectrum be- 
tween the template and iodine epochs. All velocities are 
measured relative to the zero-point denned by the template 
observation. Our periodogram analysis in Fig. |5| reveals a 
single broad peak around 1000 days at about the 0.05% false 
alarm level and well above the indicated 1% false alarm rate. 



Table 3. Relative radial velocities (RV) are given for HD 20782. 
Julian Dates (JD) are heliocentric. RVs are barycentric but have 
an arbitrary zero-point determined by the radial velocity of the 
template. 



JD 


RV 


error 




(-2451000) 


(m s" 1 ) 


(m s -1 ) 




35.3195 


24.3 


4.6 




236.9306 


4.8 


6.0 




527.0173 


9.2 


6.8 




630.8824 


31.7 


5.3 




768.3089 


-3.9 


5.1 




828.1107 


-4.8 


6.0 




829.2745 


-0.5 


8.3 




856.1353 


-2.9 


7.5 




919.0066 


6.5 


6.1 




919.9963 


6.3 


5.7 




983.8901 


11.8 


6.6 




1092.3044 


20.3 


4.7 




1127.2681 


19.6 


5.6 




1152.1631 


21.4 


5.1 




1187.1597 


25.5 


1.9 




1511.2066 


-0.3 


4.5 




1592.0484 


19.1 


4.5 




1654.9603 


22.3 


4.4 




1859.3054 


-197.6 


3.8 




1946.1383 


-14.9 


4.0 




1947.1225 


-12.3 


3.4 




2004.0015 


3.1 


3.6 




2044.0237 


7.6 


4.3 




2045.9607 


4.8 


3.8 




2217.2881 


12.1 


3.3 




2282.2204 


26.5 


3.8 




2398.9697 


27.8 


2.6 




2403.9607 


30.3 


4.8 




2576.3073 


-5.9 


3.0 




2632.2813 


-3.8 


3.2 




Table 4. Orbital parameters for the ■ 


companions 


to HD 18708 


and HD 20782. 








Parameter 


HD 187085 


HD 20782 


Orbital period P (d) 




986 


585.86±0.03 


Velocity amplitude K (m s~ 




17 


115±12 


Eccentricity e 




0.47 


0.92±0.03 


w (deg) 




94 


147±3 


Pcriastron Time (JD) 


2450912 2451687.1±2.5 


Msini (M Jup ) 




0.75 


1.80±0.23 


a (au) 




2.05 


1.36±0.12 


RMS (m s- 1 ) 




7.1 


5.0 



Our best fit solution yields an orbital period of 986 d, a 
semi-amplitude of stellar reflex motion (which we shorten to 
velocity amplitude) of 17 m s _1 , and an eccentricity of 0.47. 
The fit is shown in the upper part of Fig. [3] The RMS to the 
Keplerian fit is 7.07 m s _1 , yielding a xi °f 1-82. This solu- 
tion implies a minimum (M sin i) mass of 0.75 Mj up , and a 
semi-major axis of 2.05 au. However, the Keplerian curve of 
this eccentric orbit includes both a peak and a sharp drop 
that are not sampled by our velocities. This raises the pos- 
sibility that an orbit of lower eccentricity might fit the data 
nearly as well. Orbits with lower eccentricity deserve con- 
sideration because of the lower effective degrees of freedom 
harboured by such orbits. We tried other orbital solutions 
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10 100 1000 
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Figure 2. A periodogram of the 40 epochs of HD 187085 ob- 
tained at the AAT between 1998 November and 2005 October. 
The false alarm probabilities shown are based on the assumption 
of normally distributed errors. The dashed line is the 1% false 
alarm level. 

with fixed eccentricities. We find a wide range of eccentric- 
ities fit the data. For example, a fixed eccentricity of 0.1 
yields the same value of RMS (shown in the lower part of 
Fig. [3] This is consistent with the numerical simulations re- 
ported in Butler et al. (2000) which indicate the uncertainty 
in determining eccentricity becomes large for low amplitude 
signals. The next periastron passage for HD 187085 is 2007 
June at which time the eccentricity solution can be better 
constrained. 

Based on associating an object with others that reside 
near it on an H-R diagram and finding the best-fit rela- 
tionship between observed jitter and activity, Wright (2005) 
predict a jitter of 2.8 m s" 1 for HD 187085. Combining 
the expected error of 3 m s' 1 in quadrature with the ex- 
pected jitter error the expected RMS error of 4.25 m s _1 
is rather lower than our plausible orbital solutions which 
have ~7 m s" 1 RMS. We expect these higher RMS val- 
ues arise because the stellar jitter relationship for inactive 
stars is relatively poorly determined and also because our 
precision for HD 187085 is limited by the slightly smaller 
equivalent widths of its stellar lines relative to later spectral 
types. Thus with our observing exposures adjusted to give 
S/N=200 per exposure our precision is probably limited to 
around 4 rather than 3 ms" 1 . The lack of any observed vari- 
ation in chromospheric activity between Henry et al. (1996) 
and Tinney et al. (2002) or significant photometric variabil- 
ity within the Hipparcos dataset gives us confidence that the 
Keplerian curve arises from an exoplanet rather than from 
long-period starspots or chromospherically active regions. 



4 ORBITAL SOLUTION FOR HD 20782 

The 29 epochs of Doppler measurements for HD 20782, ob- 
tained between 1998 August and 2005 September, are shown 
graphically in Fig. 2] an d listed in Table [3] For cases of ex- 
treme eccentricity, like HD 20782, periodograms break down. 
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Figure 3. Doppler velocities obtained for HD 187085 from 1998 
August to 2005 October. The top plot shows the best fit eccen- 
tricity, the bottom plot shows a similar quality fit for a fixed 
eccentricity of 0.1. In both cases the solid lines indicate the best 
fit Keplerian with the parameters shown. 



The best-fit orbit is found by searching over a wide-swath 
of period space. In the case of HD 20782 this is simplified 
by the large velocity amplitude of the velocity variations. 

The data are well-fit by a Keplerian curve which yields 
an orbital period of 585.86 ± 0.03 d, a velocity amplitude of 
115 ± 12 m s _1 and an eccentricity of 0.92 ± 0.03. The min- 
imum (Msini ) mass of the planet is 1.77±0.22 Mj up and 
the semi-major axis is 2.26 ± 0.46 au. The RMS to the Kep- 
lerian fit is 5.0 m s _1 , yielding a reduced xl of 1.01. Wright 
(2005) predicts a jitter of 4.1 m s" 1 . Our errors come from 
Monte-Carlo simulations where the best fit Keplerian orbit, 
is recomputed based on the noise in the data. Combining our 
expected error of 3 m s" 1 in quadrature with the expected 
jitter gives an expected RMS of 5.1 m s _1 , consistent with 
the measured value. Furthermore, the lack of any observed 
chromospheric activity or photometric variations gives us 
confidence that the Keplerian curve arises from an exoplanet 
rather than from long-period starspots or chromospherically 
active regions. The properties of the candidate exoplanet in 
orbit around HD 20782 are summarised in Table [I] 

While we are confident about our procedures and the 
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Figure 4. Doppler velocities obtained for HD 20782 from 1998 
August to 2005 September. The solid line is a best fit Keplerian 
with the parameters shown in Table[I] The RMS of the velocities 
about the fit is 5.0 m s — 1 . Based on a primary mass of 0.95 Mq, 
the minimum (M sin i) mass of the companion is 1.77 Mjtjp and 
the semimajor axis is 1.36 au. 



consistency of our solution, we find that within our errors 
HD 20782 shares the status of 'most highly eccentric exo- 
planet' along with HD 80606 (e = 0.9330±0.0067, Butler 
et al. 2006; e = 0.927±0.012, Naef et al. 2001). Thus it is 
important that we examine the reality and our reliance on 
our single periastron data point (JD - 2451000 = 1859.3054, 
RV = -197.6m s _1 ). If the time of this epoch were recorded 
wrong by about 3 hours, the barycentric correction would 
be off by about 200m s _1 due to diurnal rotation. A check 
of the continuity of the logsheets rules this out and reveals 
that the periastron observations were taken under a clear 
sky. However, if the star were by chance observed near the 
Moon or in twilight the G2V reflected solar spectrum might 
drag the spectral line centroid off. 

At the time of observation the Moon as seen from Siding 
Spring was at 17 h 03 m 32 B -23° 59 57" (J2000 - from the JPL 
Horizons System, http://ssd.jpl.nasa.gov/horizons.html). 
HD 20782 is at 03 h 20 m 04 s -28°5l'l5" (J2000). Thus 
HD 20782 was nowhere near the Moon. In fact, the star is 
well off the ecliptic so can never get very close to the Moon. 
The observation time (the mid-point of the exposure) was 
four minutes before the start of morning -18° (astronomi- 
cal) twilight. At this time there should be no significant twi- 
light contamination in the two ten minute exposures taken 
of HD 20782, indeed AAPS observations routinely continue 
until -12° twilight. We are reassured of our strategy because 
based on our test exposures with UCLES looking at blank 
sky 20 minutes after -12° twilight no sky photons are de- 
tected in a 1 min exposure (above read noise) . Any solar con- 
tamination would generate a high xt statistic for the quality 
of the doppler fit. The formal internal errors for the two ob- 
servations, made consecutively at the periastron epoch, were 
5.1 and 5.5 m/s, compared with a median of 5.3 m/s for all 
observations of HD 20782. So, the individual spectra mak- 
ing up our periastron data point show no extra scatter in 
velocity. 

Despite our confidence in these two data points com- 
prising the JD - 2451000 = 1859.3054 epoch we also recom- 



puted the Keplerian fit after removing them from the data 
set. The resulting best-fit parameters are period = 599.82 d, 
periastron time (JD) = 51067.098, eccentricity = 0.732, 
to — 119.6 deg, velocity amplitude = 32.7 m s -1 and plan- 
etary mass = 0.92 Afj up . With this new fit the period is the 
same as before, but the values of eccentricity and velocity 
amplitude (and thus planetary mass) are significantly lower. 
So, although the qualitative sense of the orbit remains nearly 
the same, the epoch of low velocity data do force the veloc- 
ity amplitude, eccentricity and planetary mass to be smaller 
than if we didn't have them. Thus, while we are confident of 
the extreme eccentricity of HD 20782b, we would like to ob- 
tain further data points near periastron. At the last perias- 
tron (2005 March 27) we were unable to observe HD 20782. 
The constrained blocks of time awarded to the AAPS re- 
stricts our ability to reactively schedule time-critical obser- 
vations such as those we wish to make for HD 20782. 



5 DISCUSSION 

The two exoplanets presented here have relatively large ec- 
centricities and small velocity amplitudes. With an eccen- 
tricity of 0.925±0.030, HD 20782 has a comparable eccen- 
tricity to HD 80606 (e = 0.9330±0.0067, Butler et al. 2006; 
e = 0.927T0.012, Naef et al. 2001) which has the highest 
value of eccentricity reported to date. The companion to 
HD 187085 has a more modest and uncertain eccentricity. 
Nontheless its combination of large eccentricity and small 
velocity amplitude make it noteworthy. It may have one of 
the highest known eccentricities, alternatively, low eccentric- 
ity fits yield velocity amplitude solutions less than 15 m s _1 
making it one of only a handful of low velocity amplitude, 
long period exoplanet detections. On the other hand, if a low 
value of eccentricity is appropriate for the orbital solution 
of HD 187085 then velocity amplitude solutions of less than 
15 m s _1 will be fit placing the companion to HD 187085 
with only a handful of low velocity amplitude, long period 
exoplanet detections. 

The relative rarity of such high eccentricity planets can 
be seen in Fig. [5] The plot shows a clear deficit of eccen- 
tric planets with small semimajor axes. For planets with 
periastron distances less than 0.1 au this deficit is expected 
(e.g., Rasio & Ford 1996) and seen to be due to the tidal 
circularisation of close-in planets (e.g. Halbwachs, Mayor & 
Udry 2005). However, for periastron distances of 0.1 au the 
tidal circularisation timescale is already likely to be many 
Gyr, for inferred Q values (based on planets having semi- 
major axes with less than 0.1 au, where Q is the specific 
dissipation function). The circularisation timescale, i C i rc , is 
a sensitive function of the radius of the planet, i? p i an et (icirc 
oc ^pkLeti Goldreich & Soter 1966) and supposedly occurs 
due to dissipation in the planet, not the star (e.g., Gu et al. 
2004). Perhaps, the larger planet radii during contraction, 
in the first 10 Myr or so, helps to shorten the circularisation 
time. Thus the rising mean eccentricities versus semi-major 
axis apparent in Fig.|S]from 0.1 to 0.3 au may be explained 
by the decreasing effectiveness of the tidal circularisation 
mechanism. 

Fig. |S] in particular the lower binned plot, indicates 
that mean eccentricities rise steeply reaching a plateau be- 
yond around 0.3 au. Thereafter, the eccentricity of an ex- 
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Figure 5. The upper plot shows a scatter plot for eccentricities 
against semimajor axis for exoplanets, within 200 pc, with semi- 
major axes <10 au and masses less than 20 Mj up as given by 
table 2 of Butler et al. (2006). The values for the candidate ex- 
oplanets presented here are shown as triangles. The lower plot 
shows mean eccentricities plotted against semimajor axis for exo- 
planets. Tidal circularisation dominates the close-in planets with 
eccentricities reaching a plateau by around 0.2 au. The error bars 
are from yj (number) statistics and are only indicative. The plots 
show results from many surveys and are thus drawn from an in- 
homogeneous sample. 



oplanet appears roughly constant. Beyond 2 au there may 
be a slight deficit of high eccentricity exoplanets apparent 
from the lack of objects in the top right of the upper part 
of Fig. In order to detect highly eccentric exoplanets with 
long periods, it is necessary to have good time sampling. 
For example, the detection of the ultra-sharp periastron pas- 
sages of eccentric orbits requires relatively large numbers of 
observations as well as fortitous sampling. In the case of 
HD 20782, we only have one epoch near periastron passage. 
Nonetheless, our robust observing and data reduction pro- 
cedures and the relatively large (115 m/s) amplitude signal 
gives us confidence that we have detected a high eccentricity 
orbit. Cumming (2004) illustrates the difficulties of detect- 
ing planets of high eccentricities with limited sampling. His 
work shows that even with high signal-to-noise and repre- 
sentative numbers of data points, "there remain significant 
selection effects against eccentric orbits for e > 0.6." 

While time sampling is crucial for adequate sampling 
of highly eccentric orbits, it is also appropriate to comment 
on precision. The top left of Fig.|S|is relatively poorly popu- 
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Figure 6. The upper plot shows a scatter plot for eccentricities 
against velocity amplitude for exoplanets, within 200 pc, with 
semimajor axes 0.5—10 au and masses less than 20 Afj up as given 
by table 2 of Butler et al. (2006). The values for the exoplanets. 
The lower plot shows mean eccentricities plotted against velocity 
amplitudes. The error bars are yj (number) statistics and are only 
indicative. The plots show results from many surveys and are thus 
drawn from an inhomogeneous sample. 



lated. As discussed above this is likely to be due to difficul- 
ties associated with the detection of high eccentricity orbits. 
In addition, there maybe a relatively smaller population of 
high eccentricity exoplanets.. However, at more modest ec- 
centrcities, there is an apparent deficit of low velocity ampli- 
tude exoplanets. It is interesting to look at the objects that 
define the upper envelope of the top left of the upper part of 
Fig. |S| These objects represent the highest eccentricity ob- 
jects that have been found for a given velocity amplitude. In 
addition to the companion to HD 20782 reported here, the 
other two objects are HD 45350 b (Marcy et al. 2005a) and 
HD 11964 b (Butler et al. 2006), both recently discovered us- 
ing Keck data. Both the AAPS and Keck planet search oper- 
ate with demonstrated long-term precisions of around 3 m/s 
(e.g. fig. 1, Tinney et al. 2005 and fig, 2, Vogt et al. 2005). 
Such long-term precisions may now be rivalled over shorter 
baselines by the HARPS (e.g., Lovis et al. 2005) and Hobby- 
Eberly (e.g., McArthur et al. 2004) planet search projects 
though many of the objects plotted in Fig.l^were discovered 
by searches with lesser long term precision. While the impor- 
tance of precision is widely recognised its specific importance 
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Figure 7. The upper plot shows a scatter plot for eccentricities 
against Afsini Mj up for exoplanets within 200 pc with semimajor 
axes 0.5—10 au and masses less than 20 Mj up as given by table 2 of 
Butler et al. (2006). The values for the exoplanets presented here 
are shown as triangles. The error bars are yj (number) statistics 
and are only indicative. The plots show results from many surveys 
and are thus drawn from an inhomogeneous sample. 



for the detection of high eccentricities has been anticipated 
in the simulations of Cumming (2004) and Halbwachs et 
al. (2005). Cumming's fig. 4 shows that the detectability of 
exoplanets with e > 0.7 is a function of precision. Thus, 
simulations and the exoplanets detected so far serve to sug- 
gest that radial velocities with good long-term precision are 
important to robustly sample the parameter space occupied 
by highly eccentric planets. 

The observational difficulties inherent in the detection 
of low velocity amplitude - high eccentricity planets pre- 
dicted by Cumming and suggested in the upper part of Fig. El 
appear to be borne out by the binned (lower part) of the fig- 
ure. These binned data suggest that the velocity amplitude 
of exoplanets increases with mean eccentricity in Fig. al- 
beit with low confidence since it is almost possible to draw a 
flat line through all the error bars. While we have attempted 
to remove the effects of circularisation processes from Fig. 
|S| (by removing all planets with semi-major axes less than 
0.5 au), it is plausible that the suggested trend may have 
a physical origin. In particular, it is notable that beyond 
the tidal circularisation limit (usually assumed to be 0.1 au) 
spectroscopic binaries have significantly higher eccentricities 



than exoplanets (Halbwachs et al. 2005) and that high mass 
planets (also having high velocity amplitudes) have system- 
atically higher eccentricities than low mass planets, e.g., 70 
Vir (Marcy & Butler 1996) and fig. 5, Marcy et al. (2005b). 
Fig.[7|shows this correlation of eccentricity with mass for the 
exoplanet sample of Butler et al. (2006). For large velocity 
amplitudes and large mass planets this cannot be a selec- 
tion effect nor can it be caused by errors because the most 
massive planets induce the largest Keplerian amplitudes, al- 
lowing precise determination of eccentricity. To quantify the 
reality of the suggested trends, to assess observational bias 
at high eccentricities and to recover the underlying eccentric- 
ity distribution of exoplanets analysis of the radial velocity 
data for each different survey is required. Such work has been 
started by Cumming et al. (1999, 2003) and independently 
by the Anglo- Australian Planet Search. 

By the standards of our Solar System, the majority of 
exoplanets beyond 1 au can be seen from Fig.|S]to be in ec- 
centric orbits. This contrast with the Solar System may be 
at least partly explained by as yet unseen long-period com- 
panions. Indeed the fits for a few of the long-period eccentric 
planets do already include longer-period trends though these 
are generally rather modest and sometimes uncertain, e.g., 
our 'best fit' solution for HD 187085 includes an additional 
slope of 1.3±1.0 m s" 1 . So while individual cases require 
careful consideration, the lack of substantial trends in the 
current sample of high eccentricty exoplanets (e.g., Butler 
et al. 2006) means that the majority of the observed eccen- 
tricities need to be explained by high formation eccentricties 
or by appropriate eccentricity pumping mechanisms. Since 
the derived masses of exoplanets are directly proportional 
to measured velocity amplitude it is important to account 
for the likely observational eccentricity bias when investigat- 
ing mass dependencies of exoplanets and in particular when 
trying to understand the relative importance of different ec- 
centricity pumping mechanisms (e.g., Takeda & Rasio 2005). 
We are encouraged to proceed in earnest with our project to 
examine the detectability of Anglo- Australian Planet Search 
exoplanets of different eccentricities (and other parameters), 
given our errors and time sampling. This should improve our 
search and will put us in a better position to better under- 
stand the underlying distribution of eccentricity with mass 
and semimajor axis for exoplanets. 
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